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The first direct measurement of 12 C( 12 C,n) 23 Mg at stellar energies 
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Neutrons produced by the carbon fusion reaction 12 C( 12 C,n) 23 Mg play an important role in stellar 
nucleosynthesis. However, past studies have shown large discrepancies between experimental data 
and theory, leading to an uncertain cross section extrapolation at astrophysical energies. We present 
the first direct measurement that extends deep into the astrophysical energy range along with a 
new and improved extrapolation technique based on experimental data from the mirror reaction 
12 C( 12 C,p) 23 Na. The new reaction rate has been determined with a well-defined uncertainty that 
exceeds the precision required by astrophysics models. Using our constrained rate, we find that 
12 C( 12 C,n) 23 Mg is crucial to the production of Na and A1 in Pop-Ill Pair Instability Supernovae. 

It also plays a non-negligible role in the production of weak s-process elements as well as in the 
production of the important galactic q-emitter eo Fe. 

PACS numbers: 26.20.-f, 26.20.Np, 26.30.-k, 25.70.Jj, 26.20.Kn 


The first stars in the early Universe formed about 400 
million years after the big bang. Verification of the ex¬ 
istence of these stars is important for our understand¬ 
ing of the evolution of the Universe [I]. It has been 
predicted that for Population-Ill (metal-free stars [2j) 
stellar production yields, the abundances of odd-Z el¬ 
ements are remarkably deficient compared to their ad¬ 
jacent even-Z elements fj. Astronomers are searching 
for long-lived, low mass stars with the unique nucleosyn- 
thetic pattern matching the predicted yields j3j. The 
relevance of 12 C( 12 C,n) 23 Mg in the first stars has been 
discussed by Woosley, Heger, and Weaver [5]. By the 
end of helium burning in Pop-Ill stars, the neutron to 
proton ratio in the ash is almost exactly 1. However, in 
the subsequent carbon burning phase, frequent /? + decay 
of produced 23 Mg converts protons into neutrons, thus 
increasing the neutron to proton ratio. A slight excess 
of neutrons would significantly affect the abundances of 
the odd-Z isotopes with neutron to proton ratios higher 
than 1, e.g. 23 Na and 27 Al. 

12 C( 12 C,n) 23 Mg is also a potentially important neu¬ 
tron source for the so-called weak s-process occurring in 
massive Pop-I (metal-rich |3] ) and Pop-II (metal-poor [5] ) 
stars. The weak s-process takes place during the core 
helium and shell carbon burning phases and is largely 
responsible for the s-process abundances up to A«90 [6]. 
Pignatari et al. recently performed a study of the weak 


s-process during carbon shell burning for a 25 Mq stel¬ 
lar model using different 12 C( 12 C,n) 23 Mg rates [7 . They 
found that a factor of 2 precision or better would be de¬ 
sirable to limit its impact on the s-process predictions to 
within 10%. 

Stellar carbon burning has three main reaction chan¬ 
nels: 

12 C + 12 C -A 23 Mg + n - 2.60 MeV 
-A 23 Na + p -(- 2.24 MeV 
-A 20 Ne + a + 4.62 MeV 

With Q < 0, the probability of decay through the neutron 
channel is weakest among the three at the low energies 
relevant for astrophysics. For a typical carbon shell burn¬ 
ing temperature Tg = 1.1, the important energy range 
for this channel is 2.7 < E cm < 3.6 MeV. The reaction was 
first studied in 1969 by Patterson et al. jB] who measured 
the cross section over the range E cm =4.23 to 8.74 MeV 
by counting /3-rays from 23 Mg decays. From this mea¬ 
surement, a constant neutron branching ratio, /3 n =2%, 
was deduced [§]. Later Dayras et al. extended the mea¬ 
surement down to E cm =3.54 MeV by counting the q-rays 
emitted following the 23 Mg beta decay. The experimen¬ 
tal uncertainty is about 40% at E cm ss 3.8 MeV and in¬ 
creases to 90% at the lowest energy [10 . To estimate the 
cross section at the stellar burning energies, Dayras et al. 
had to rely on an extrapolation of the experimental data 
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based on a Hauser-Feshbach statistical model calculation 
m Because of the unique molecular resonances exist¬ 
ing in the 12 C+ 12 C fusion reaction pj2], their calculation 
could only be renormalized to the average trend of the 
data while the resonant behavior of the 12 C+ 12 C fusion 
reaction was ignored. The maximum deviation between 
the experimental result and the renormalized statistical 
model prediction is more than a factor of 4 (see Fig. [lj. 
Nevertheless, based on the statistical model extrapola¬ 
tion, this work recommended a neutron branching ratio 
of /?„= 0.011%, 0.11%, 0.40% and 5.4% at T 9 =0.8, 1.0, 
1.2, and 5, respectively, though no attempt was made to 
quantify the uncertainties in these predictions jT0] . 

In 1988, Caughlan and Fowler (CF88) excluded this 
result from their rate compilations m ■ Instead, they 
recommended /3„= 0 (T 9 < 1.75), /3„= 5% (1.75 < T 9 < 3.3) 
and (3 n = 7% (3.3 < Xg < 6.0). This rate was adopted by 
REACLIB after fitting the CF88 ratio with the standard 
REACLIB formula [14]. Pignatari et al. attempted to 
use the Dayras rate |7j, however it was later discovered 
that the analytic formula for /3 n taken from the paper m 
contained a typographical error resulting in a significant 
deviation from the intended value below Tg =1.5 |15| . So 
far, to our knowledge, the correct Dayras rate has only 
been implemented in the stellar code KEPLER [m nz]. 

In the following, we report on the first direct measure¬ 
ment of this reaction into the stellar energy range as well 
as an improved method for extrapolating the experimen¬ 
tal results through the remaining unmeasured energies 
relevant for carbon shell burning. Based on the new ex¬ 
perimental result, a new reaction rate is recommended 
together with a well-defined uncertainty. The impact on 
the nucleosynthesis in massive stars is also discussed. 

The experimental work was performed at the Univer¬ 
sity of Notre Dame’s Nuclear Science Laboratory using 
the 11 MV FN tandem Van de Graaff accelerator. Car¬ 
bon beams were produced at energies ranging from 5.1 
to 8.7 McV (lab frame) with typical currents on tar¬ 
get between 0.5 and 1.5 p/iA. The beam energy calibra¬ 
tion was checked by measuring the reaction thresholds 
of 7 Li(p,n) and 19 F(p,n) as well as 12 C(p,p) resonant 
scattering |18] . The maximum energy deviation was less 
than 0.1%. A 1-mm thick hydrogen-free Highly Ordered 
Pyrolytic Graphite (HOPG) target made from natural 
carbon was used to control the hydrogen-induced back¬ 
ground [19]. The target was cooled by circulating deion¬ 
ized water through the supporting flange, which was cen¬ 
tered in a block of polyethylene containing 20 3 He pro¬ 
portional counters arranged around the beam axis in two 
concentric rings mum. 

The main sources of beam-induced neutron back¬ 
ground were from the reactions 13 C( 12 C,n) 24 Mg and, 
to a lesser extent, 2 H( 12 C, 13 N)n [TH] , With a large 
positive Q-value (8.99 MeV) and the relatively high 
natural abundance of 13 C in the target (1.1%), neu¬ 
trons from 13 C( 12 C,n) 24 Mg dominate the total yield at 



(MeV) 


FIG. 1. (Color online) Upper: The 12 C( 12 C,n) 23 Mg S*-factor 
results from the present measurement (black squares) com¬ 
pared with previous data sets from Dayras 1977 [TD] (purple 
circles) and Bucher 2013 [151 l2ll (magenta triangles). Also 
shown is the neutron branching ratio calculated by Dayras 
applied to the total 12 C+ 12 C fusion S*-factor recommended 
by CF88 p[3] (red solid line) and the new extrapolation from 
this work (blue circles). Only statistical errors are shown for 
the experimental data, while the extrapolation includes both 
statistical and a 40% systematic error. Lower: The integrand 
from Eq. fusing the various data sets is plotted. 

very low beam energies approaching the 12 C( 12 C,n) 23 Mg 
reaction threshold. To estimate its contribution, the 
13 C( 12 C,n) 24 Mg reaction was studied with the same 
setup using a 13 C beam with energies ranging between 
9.5 and 5.4 MeV. Since the cross section for this reac¬ 
tion is much higher, relatively low beam intensities (~50 
pnA) with shorter run times were sufficient. The nor¬ 
malized 12 C( 13 C,n) 24 Mg yield was then subtracted from 
the measured total neutron yield recorded with the 12 C 
beam m- 

The room background rate was measured to be 
9.015(92) evts/min, which dominated the yield at ener¬ 
gies below E cm =3.0 MeV. The background contribution 
from 2 H( 12 C, 13 N)n was studied using a thin TiDg target 
with thick Cu backing. After removing the room back¬ 
ground, this contribution was found to be less than 5% 
of the total yield at E cm =3.3 MeV increasing to 19% at 
3.1 MeV. 

The detector efficiency has been simulated using 
Geant4 and MCNP in the range 0 to 3 MeV and experi¬ 
mentally validated with monoenergetic neutrons between 
50 and 650 keV [2D]. For this experiment, the Geant4 
code was modified to include the strong angular depen¬ 
dence of the neutron energy from the 12 C( 12 C,n) 23 Mg 
kinematics. An isotropic angular distribution in the 
center-of-mass frame was assumed for the neutrons pro¬ 
duced by 12 C( 12 C,n) 23 Mg and 13 C( 12 C,n) 24 Mg. The ef¬ 
ficiency was found to vary smoothly between 30% and 
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50% as the average neutron energy decreases with beam 
energy [2Tj . To check the effect of the assumed angular 
distribution, we changed the isotropic angular distribu¬ 
tion to match 2 nd and 4 th order Legendre polynomials 
122) . and found a nearly constant relative drop in effi¬ 
ciency of 9% and 5%, respectively, in the range of E cm = 3 
to 5 MeV. Since our experiment does not measure angu¬ 
lar distribution, a ±5% systematic uncertainty has been 
assigned for this effect. To test our efficiency calculation, 
we measured the 12 C( 12 C,n) 23 Mg cross section indepen¬ 
dently by detecting the activity of 23 Mg [15], ;2Tj. The 
results gave good agreement in the overlapping energy 
range (as shown in Fig. [l]). 

The cross section for the 12 C( 12 C,n) 23 Mg reaction was 
determined by differentiating the thick target yield [23] , 
In Fig. [l] it has been converted to a modified S-factor 
(S*) for comparison with previous results. S* differs from 
the standard S-factor only by the multiplier, exp(0.46 E), 
which is added to account for the finite size of the nucleus 
[5] . It is seen that the new results display good agreement 
with previous measurements in the overlapping energy re¬ 
gion while extending much deeper into the astrophysical 
energy range. A new resonance at E cm = 3.4 MeV is ob¬ 
served in the neutron channel. This resonance was also 
observed in earlier measurements of the proton and alpha 
channels mciMi- Our measurement includes a 15% 
systematic uncertainty which primarily results from the 
uncertainties in the beam current (10%), beam energy 
(2%), detector efficiency (6%) [20], angular distribution 
(5%), and stopping power (7%) (21, 28]. The Dayras re¬ 
sults also have an additional systematic uncertainty of 
16% [TO] not shown in Fig. [TJ 

An extrapolation is required to estimate the reaction 
cross section at the lower energies beyond experimental 
reach. As mentioned earlier, Dayras et al. provided a 
renormalized statistical model calculation for this pur¬ 
pose. However, the large discrepancy between the ex¬ 
perimental data and their theory calls into question the 
reliability of the extrapolation. To provide a better pre¬ 
diction including the effect of the molecular resonances in 
the entrance channel, a novel extrapolation method has 
been developed based on experimental information from 
the mirror reaction 12 C( 12 C,p) 23 Na. The predicted neu¬ 
tron cross section, cx n ( pre d), is obtained using the formula 

N 

. _ V" a n, (th) _ 

^n(pred) — / ^pj(exp) 

7^0 CT P*(th) 

where N is the highest available decay channel in the 
residual 23 Mg, which depends on the reaction energy. For 
E cm <4.6 MeV, only the no and ni channels are open. 
The theoretical ratio, <r ni (th) /<r Pi (th), is calculated using 
TALYS [29] combined with entrance channel spin popu¬ 
lations supplied from a coupled-channels calculation by 
Esbensen j3U]. The resonances in 12 C( 12 C,ni) 23 Mg and 
12 C( 12 C,Pi) 23 Na originate from both the molecular res¬ 


onances in the entrance channel and the characteristic 
resonances in the final decay channels. The traditional 
statistical model calculation employed by Dayras uses the 
optical model and assumes a high level density to describe 
the entrance and exit channels and therefore could only 
reproduce the average trend of the experimental data. 
In our approach, the complicated molecular resonance 
associated with the entrance channel is embodied in the 
experimental cross sections (cr p .( exp )) of 12 C( 12 C,Pi) 23 Na, 
the mirror system of 12 C( 12 C,n i ) 23 Mg, while the statis¬ 
tical model is only used to predict the decay width ratio 
between the n^ and p^ channels. Since the proton energy 
resolution in the Zickefoose experiment from Ref. [26] 
was insufficient to resolve po from pi, only the sum, 
<7 po + <r pi , is available for E cm <4 MeV. Eq. [I] has been 
modified to accommodate the combination of po and pi. 
Additionally, the measurements of 12 C( 12 C,pi) 23 Na by 
Fang et al. [27] performed at Notre Dame in the energy 
range 3<E cm <6 MeV have also been used to predict 
the 12 C( 12 C,n) 23 Mg cross section [T3]. In this case, up 
to N=6 possible decay channels are required for the pre¬ 
diction calculated in Eq. [l] 

Figure [2] shows the ratios between our measured 
12 C( 12 C,n) 23 Mg cross section cr„( exp ) and the two cr„( pre d) 
based on the Zickefoose and Fang proton data sets plot¬ 
ted as a function of E cm . The average ratios (standard 
deviations) for the Zickefoose- and Fang-based predic¬ 
tions are 0.9(4) and 0.9(3), respectively. The ratios to 
the Dayras calculation are also shown for comparison. 
The large deviation at E cm ~4.8 MeV has been elim¬ 
inated by our approach. The fluctuations, which are 
larger than the quoted statistical uncertainties, reflect 
the systematic errors associated with our extrapolation. 
They consist of the systematic errors in the proton mea¬ 
surements, the assumed entrance channel spin popula¬ 
tions, and the TALYS calculation used in the prediction 
of cr ni (th)/t J Pi (th)- To provide better consistency with 
the experimental 12 C( 12 C,n) 23 Mg data, our extrapola¬ 
tion has been renormalized by the factor 0.9. We have 
adopted 0.4 as the systematic error in accordance with 
the Zickefoose-based prediction since that data set was 
used for the extrapolation, being the only one to reach 
sufficiently low energies. 

The new cross section defined by our extrapolation 
and experimental data has been used to calculate the 
12 C( 12 C,n) 23 Mg reaction rate by the following equation: 

/ 8 \ V2 1 r 00 / E \ 

m =y ^ ) iE < 2 > 

To highlight the important stellar energy range for a 
typical carbon shell burning temperature Xg= 1.1, the in¬ 
tegrand of Eq. [2] (Gamow yield) is computed and shown 
in Fig. [T| Our measurement covers about half of the stel¬ 
lar energy range. It reduces the dependence on extrap¬ 
olation in the astrophysical reaction rate and provides a 
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FIG. 2. (Color online) The ratio of our 12 C( 12 C,n) 23 Mg 
cross section data a n ( exp ) to the two <7 n (pred) based on the 
12 C( 12 C,p) 23 Na data from Zickefoose [26] (blue diamonds) 
and from Fang et al. [ 27 j (red squares). The solid black line 
shows the average ratio of the Zickefoose data below 4.0 MeV 
while the dashed lines represent 1 standard deviation. As a 
comparison, the ratios of our cr n ( exp ) to the Dayras prediction 
are shown as black circles. 
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FIG. 3. (Color online) The ratios of various published 
12 C( 12 C,n) 23 Mg rates to the one determined in this work. The 
uncertainties of our new rate are indicated by the red lines. 
The rate determined by Dayras is in good agreement with our 
new rate. The erroneous Dayras rate, stemming from a typo 
in the publication, is labeled DayrasT. For comparison, we 
also show the rates from CF88 m and two separate editions 
of REACLIB (from Thielemann et al. [14] and JINA [3l'|). 


base for examining the systematic uncertainty of extrap¬ 
olation. 

Because of its endothermic character, a good fit of the 
12 C( 12 C,n) 23 Mg reaction rate was difficult to achieve us¬ 
ing the standard REACLIB format. Following the con¬ 
vention of Dayras m, our 12 C( 12 C,n) 23 Mg rate has been 
normalized to the standard CF88 12 C+ 12 C total fusion 
rate. The neutron branching ratio, fi n , has been fitted 
and listed below. 


fin 


0.11954 exp 


(T 9 < 1.5) 


0.16446 2.57495 

-o-1- n -b 

t 9 3 t 9 2 


0.2212 [1 - exp(—0.13597 T g + 0.158)] 


1.94145 \ 
T9 ) 


(1.5 < T g < 2.5) 

0.048811 [l - exp(—2.1124 T g + 3.8791)] 



Atomic Number 

FIG. 4. (Color online) The logarithmic ratio of elemental 
yields with 12 C( 12 C,n) 23 Mg to those without 12 C( 12 C,n) 23 Mg 
for 18 Mg (top) and 200 Mg (bottom) Pop-Ill stars. The red 
points and their error bars correspond to the result obtained 
with the average and upper/lower limits determined by this 
work.The result obtained with the Dayras rate is shown as 
black points. 


(2.5 <T g < 5.0) 

= 0.04875 (T 9 > 5.0) 

The uncertainty for the reaction rate is estimated 
based on the error bars of experimental and extrapolated 
cross sections. Comparisons among the existing reaction 
rates are shown in Fig. [3] The various rates have been 
plotted as a ratio to our rate in order to compare them on 
a linear scale over a large temperature range. It is seen 
that only the Dayras rate agrees with our new rate within 
the quoted uncertainty. At typical carbon shell burning 
temperatures T 9 ~l.l-1.3, the uncertainty is less than 
40% which is sufficient for studying the weak s-process. 
The uncertainty is reduced to 20% at T 9 ~ 1.9-2.1 which 
is relevant for explosive carbon burning. 

The impact of 12 C( 12 C,n) 23 Mg on the nucleosynthetic 
pattern of a 200 M 0 Pair Instability Supernovae (PI SNe) 
has been investigated using the ID stellar evolution code, 
KEPLER 1161 !?]■ The ratio of the production yields 
with and without 12 C( 12 C,n) 23 Mg is shown in Fig. [ 4 ] 
It is clear that this reaction is important for the nucle¬ 
osynthesis of odd-Z elements such as F, Na, and Al. By 
including our 12 C( 12 C,n) 23 Mg rate in the calculation, the 
production of 23 Na is increased by a factor of 5 (0.7 dex) 
with an uncertainty less than ~10%. The yield of 2, A1 
is increased by nearly a factor of 2 (0.3 dex). We have 
also explored the impact of 12 C( 12 C,n) 23 Mg on an 18 M 0 
Pop-Ill star. A moderate enhancement of up to 30% is 
found for odd-Z elements. 

The role of 12 C( 12 C,n) 23 Mg in Pop-I stars has been 
overlooked in most studies of the weak s-process because 
both CF88 and REACLIB (from Thielemann )14| ) essen¬ 
tially turn it off at carbon shell burning temperatures. 
To illustrate its impact on the nucleosynthesis in massive 
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stars, an 18 M 0 Pop-I star has been investigated using 
KEPLER with two different scenarios: including and ex¬ 
cluding 12 C( 12 C,n) 23 Mg. By comparing these two pro¬ 
duction yields, an enhancement of ~ 10 % is found for a 
number of weak s-process isotopes, such as 70 Zn, 76 Ge, 
82 Se, 86 Kr, 85,87 Rb and 96 Zr, arising from the additional 
neutron production from 12 C( 12 C,n) 23 Mg. The origins of 
these elements are rather complicated including He-, C- 
, Ne-burning, s-process in AGB stars and the r-process. 
Even within the weak s-process there are a number of 
uncertainties that can affect the final abundance pattern 
PH IMS . Our result clarifies the ambiguities associ¬ 
ated with the 12 C( 12 C,n) 23 Mg rate. Furthermore, about 
10% enhancements are also observed for 46 Ca and 60 Fe. 
The production mechanism of 46 Ca is important for the 
understanding of the 48 Ca/ 46 Ca anomaly in meteorites 
[55] . while the production of 60 Fe is an important topic 
in y-ray astronomy (36j . 

In summary, we have measured the 12 C( 12 C,n) 23 Mg 
cross section for the first time within the Garnow window 
for the stellar carbon burning processes. Our measure¬ 
ment covers half of the important energy range. For the 
lower unmeasured energies, we have developed a novel ex¬ 
trapolation method based on the 12 C( 12 C,p) 23 Na chan¬ 
nel. A new reaction rate has been determined with, for 
the first time, a quantified uncertainty that satisfies the 
precision required from astrophysics models. As a result, 
the ambiguity arising from the uncertain 12 C( 12 C,n) 23 Mg 
reaction rate has been eliminated. With our new rate, we 
find that 12 C( 12 C,n) 23 Mg is crucial for constraining the 
production of Na and A1 in Pop-Ill Pair Instability Su¬ 
pernovae, and it plays a non-negligible role in the produc¬ 
tion of weak s-process elements as well as the production 
of the important galactic 7 -emitter 60 Fe. 
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